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We transform the two-dimensional Dirac-Weyl equation, which governs the charge carriers in 
graphene, into a non-linear first-order differential equation for scattering phase shift, using the so- 
called variable phase method. This allows us to utilize the Levinson Theorem to find zero-energy 
bound states created electrostatically in realistic structures. These confined states are formed at 
critical potential strengths, which leads to us posit the use of ‘optimal traps’ to combat the chiral 
tunneling found in graphene, which could be explored experimentally with an artificial network of 
point charges held above the graphene layer. We also discuss scattering on these states and find the 
m = 0 states create a dominant peak in scattering cross-section as energy tends towards the Dirac 
point energy, suggesting a dominant contribution to resistivity. 

PACS numbers: 


I. INTRODUCTION 

The electronic properties of the two-dimensional (2-D) material graphen^i^^ are of great interest due to the 
quasi-relativistic nature of it s sp ectrum. Interesting transport effects such as chiral (Klein) tunnehng,l^l^3 vacuum 
polarization,!^ atomic collaps^^E^ and the minimum conductivity at the Dirac pointP^ have been widely discussed in 
the literature. The topic of elastic scattering in clean, low-temperatur e grap hene, which can occur due to charged 
impurities, ripples or strain fields, has been addressed by many authors 

However, despite its extraordinary properties there is a major obstacle stopping graphene from usurping silicon 
in the electronics industry, namely is its lack of a bandgap. This frustrates attempts to perform digital logic with 
graphene due to the difficulty in turning off the flow of chiral charge carriers which always wish to conduct. Attempts 
at opening a gap in monolayer graphene have focused on cutting into nanoribbons,!^ chemical functionalizatioiP^ and 
strain engineering,^^ which can unfortunately blunt the remarkable electronic properties which makes graphene so 
attractive in the first place. Here we propose a method not to open a gap, but to switch off the chiral tunneling by 
considering zero-energy states, when the Fermi energy coincidences with t he Dir ac points, such that fully confined 
states are predicted to be able to form due to the absence of pseudospin.l^SlIIl iact, these states are the most 
important factor when considering resonant scattering in graphene. 

Despite the appearance of sophisticated experimental techniques for probing resonances and the modification of 
the density of states in the continuunJ^ the search for fully-confined (square-integrable) states remains a significant 
ongoing task. Efficient manipulation of the Fermi level requires the presence of a back-gate in close proximity to the 
graphen^ which makes the numerous beautiful results stemming from the long-range behavior of the bare Coulomb 
potentiaP to be of somewhat far from experimental reality, as the presence of image charges in the gate material (or 
screening effects) make any realistic potential fall at large distances faster than l/r.!^It is important to emphasize 
that any fast-decaying potential cannot produce a bound state at nonzero energy.!^ Indeed, the asymptotic of the 
wavefunction is a Bessel function decaying asymptotically only like and so one is led to consider zero-energy 

states instead. 

Whilst lo w-energ y resonant scattering in monolayer graphene has been intensively studied by previous authors both 
theoreticalljP^SHini and experimentallySSEll the importance of fully-confined zero-energy states in realistic structures 
has not been fully appreciated until recently.!^ Thus far, only quasi-bound states, where only one wavefunction com¬ 
ponent is confined or when t he wa vefunction is non-square integ rable, have been considered for resonant scattering. 
Previously, only circu lar well j^^ * ^^ or the Coulomb potential,!2l^have been investigated, but in this paper we concen¬ 
trate on smooth,!2fil2II short-range potentials which are defined with two parameters, characterizing both strength and 
spread. 

We study confined states and resonant scattering in graphene due to either scanning probe microscopy (SPM)P^ 
tip-induced potentials or due to some charge displaced out of the plane,!^ with careful consideration of truly bound 
zero-modes. The strength parameter can arise due to, for example, the size of the charge on the SPM, whilst the 
spread parameter is linked to the distance from the SPM tip to the graphene. We investigate both the conditions 
required for a zero-energy bound state to form and the effect of such states in our study of the energy dependence of 
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scattering cross-section and resistivity contributions of resonant scatterers. We also note there is increasing interest 
in zero-modes of the Dirac equation in the condensed matter community due to the possibilities of both observing 
the elusive Majorana fermionPSl or indeed fractionally charged excitationsP^^ 

To carry out our investigations into realistic, short-ra nge (d ue to the necessity of a gate in all measurements) 
potentials we develop the variable phase method or VPMpEU which was found to be useful for tackling scattering 
problems governed by the Schrodinger equation in 2 -D,l®SSl for use with the 2-D Dirac-Weyl equation - allowing 
us to consider the charge carriers of graphene. The VPM was originally developecP^ for use with non-relativistic 
wave equations in the 1930s as a neat tool to calculate physically relevant quantities directly, rather than having 
to extract them from the wavefunction, and has recently been developed for use with the Dirac-Weyl equation in 
quasi-one-dimensional problems due to the intense interest in graphene in the condensed matter community.^ Here 
we derive a first-order equation from which we can immediately find the scattering phase shift. This is advantageo us as 
important physical properties directly follow, such as: the number of bound states (from the Levinson Theorem) pSED 
the scattering and transport cross-sections (from standard elastic scattering theory)}^ the number of states around a 
potential barrier (using the Friedel sum rule) and the energy change due to the impurity interacting with neighboring 
electrons (the Fumi theorem) .1^ We have neglected effects due to ri ppling of or dislocations in the graphene membrane. 
We also do not discuss scattering by multiple electrostatic barrier^^^l^Sl qj. ]-,y magnetic barriers,!^^^^ but our method 
can be generalized to account for the presence of vector potentials. 

The rest of this work is as follows. We outline the logic behind the VPM and give a detailed derivation of the phase 
equation for the Dirac-Weyl equation in Sec. [ITj In Sec. m we compare the results of the VPM against tight-binding 
calculations on a finite flake and then apply elastic scattering theory in conjunction with the VPM, which allows 
us to investigate both the nature of confined states in graphene and the influence of short-range scattering on some 
of its transport properties. Finally, we summarize and discuss our results in Sec. m In Appendix we show how 
quantum confinement of zero-energy states can occur for a quite general class of potential wells decaying like a power 
law, whilst in Appendix]^ we detail the key results for the simple case of scattering by a circular finite potential well. 


II. FORMALISM 

The 2-D Dirac-Weyl Hamiltonian governing the low-energy charge carriers of graphene on a Dirac cone ifl 

H = vpa-■ p + U{r) + azA, (1) 

where up « c/300 is the Fermi velocity of the Dirac particles, cr = {crz;,cry,az) are the Pauli spin matrices, U(r) is 
a central potential and we include A as the mass term for generality. A non-zero mass can arise due to chemical 
modifications or by strain engineering.^^ We move into polar coordinates (r, 0) for circular symmetry, and separate 
the variables via the following ansatz for the two-component spinor wavefunction 


d>(r, 9) = 


( XA{r) 


V *e*®XB(?’) 


m = 0, ±1, ±2,... 


( 2 ) 


where the subscripts A and B label the two sublattices of the graphene chicken wire lattice. This choice of wavefunction 
leads to two coupled first-order differential equations for the radial wavefunction components xa,b(j’) 


d m + 1 
dr r 


dr r 


XB = {k-V{r)- A)xa, 
XA = {k-V{r) + A)xb, 


(3a) 

(3b) 


with V{r) = U{r)/hvp, A = A/hvp and k = E/Hvp, where E is the eigenenergy. Re-arranging Eqs. (§ into a 
second-order differential equation for a single radial wavefunction component XAij") we obtain for the massless case 


d^ 


-lJiXA{r)+ - + 


dr"^ 


I 


1 


dV{r) 


r k — V(r) dr 


ixAr) + Uk-Vir)f - - % I = 0- W 

dr \ r k — V [r) dr r^ 


We consider potentials of the form V(r —>■ oo) = 0 such that at large distances Eq. Q reduces td^ 


dr^' 


IZ^XAir) -f - -j-XAir) + U - ^ XAir) = 0, 


dr 


m 


( 5 ) 
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which is the Bessel equation with the well-known solution amJm{kr) — braNmikr), or equivalently 

XA(r) = Am [Jni{kr) cos(i5„) - Nm{kr) sin((5„)], (6) 

where Jm(kr) and Nm(kr) are the Bessel functions of the first and second kinds respectively, and 5m = arctan(6m/am) 
is the scattering phase shift, arising from the difference in phase of the wavefunction at r —)■ oo compared to the free 
particle case. 

We now implement the VPM by treating the constants Am and 5m as functions of the radial coordinate r, such 
that 


XA{r) = Am(r) [Jm{kr) cos{5m{r)) - Nm{kr) sin((5m(r))], 


(7) 


where Am(r) is called the amplitude function and the phase function 5m is the phase shift arising from a potential 
cut-off at a distance r. To completely define these newly introduced functions Amir) and 5m we make the following 
ansatz for the first derivative of xa{t) with respect to r 


Xa W = [J'mikr) cos((5™(r)) - N’m{kr) sin((5™(r))], (8) 

where ’ denotes differentiation with respect to r. Now, setting the direct derivative of Eq. 0 equal to Eq. ([^ suggests 
the following useful condition 


A'm{r) _ , ' Jm{kr) sin(J^(r)) -b Nmikr) cos{5m{r)) 

Amir) ^ Jmikr) cos(5m(r)) - Nmikr) sin(5m(r))' 


Upon substituting Eq. Q and Eq. 
component Xsir) is 


into the lower coupled Eq. (3b I we naturally find the lower radial wavefunction 


= k^-Vir) + ^Jmikr)^ cosi5mir)) - [-N'mikr) + sin((5^(r)) . (10) 

We can now utilize the upper coupled Eq. (3a): upon substituting in Eq. ([^ and Eq. (10) and eliminating the 
amplitude function Amir) via the application of the condition Eq. we obtain the following first-order differential 
equation 


d 

dr 


5mir) 


7rr r 1 dVir) 
[k-Vir) dr 


— p(r)^ + (U (r)^ - 2kV (r)) p(r) , 


Pir) = Jmikr) cos((5m(r)) - Nmikr) sin((5„(r)), 
= J'mikr) cos((5„(r)) - N^kr) sin((5„(r)), 


( 11 ) 


where we have introduced the auxiliary functions p(r) and qir) and in addition have used the Wronskian of the Bessel 
functions W{Jmix), Nmix)} = Jmix)Nmix) — Nmix)Jmix) = —2/(7ra;) to simplify the final expression.!^ Eq. (11) is 
the so-called phase equation, and is subject to the initial condition 5mi0) = 0, as follows from being in the free particle 
limit. We can see from Eq. 0) how the potential V (r) gradually accumulates the desired phase shift starting from 
5m iO) = 0 and finishing with the total phase shift of the scattering problem, given by 


5m = lim 5m ir). (12) 

r—>-oo 

This condition ensures the phase shift is uniquely defined, avoiding an ambiguity of tt that appears in other methods.!^ 
When investigating bound states in the massless case we can only consider zero-energy states, where the Neumann 
function is divergent and so the following condition is implied 


5m=nn, n= 1,2,3... (13) 

Eq. ( [T^ is related to the Levinson’s theorem for massless 2D Dirac particles, which states a relation between the 
phase shift at zero-momentum and the number of bound states. 

Please note when considering the massive 2D Dirac particles, e.g. as found in h-BN or gapped graphene,!^ which 
allow bound states at finite energy, an equation analogous to the phase equation Eq. 0 can be derived for treating 
confined states 


dr 


Vmir) 


-rfir) 


1 dV ir) 
k — Vir) dr 



^/(O) + {Vir)^ - 2/cU(r)) /(r) 


fir) = ImiKr) cosipmir)) - KmiKr)sinir]mir)), 
gir) = I'mi^r) cosipmir)) - iX(„(Kr) sin(??^(r)). 


(14) 
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where and Km{nr) are the modified Bessel and Neumann functions respectively and the effective wavevec- 

tor K = (A^ — Eq. (14) has been simplifiecP^ by noting the Wronskian of the modified Bessel functions 

W{Im{x)TKm{x)} = — Km{x)I'^{x) = —1/x. Notably, Eq. (14) is also relevant for considerations of the 

surface states on 3-D topological insulators such as Bi 2 Te 3 , where in this case the mass term arises from the exchange 
energy from a magnetic insulator.l^ When considering bound states we note the modified Bessel function of the first 
kind is divergent and so the following condition is implied 


= I n - - ) TT, 


n = 1, 2,3... 


(15) 


which is a representation of the Levinson theorem for the massive 2D Dirac equation.l^SIlll 


III. RESULTS 
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FIG. 1: (Color online) Plots of scaled phase shifts 5m/7r against potential strength Vod for massless Dirac particles, of energy 
tending towards zero, incident on (top) the Lorentzian potential and (bottom) the model potential of Eq. (16) with the realistic 
1/r^ decay. We show resnlts for angular momentum m = 0,1,2 corresponding to the solid line (red), dashed line (blue) and 
dotted line (green) respectively. 


We now check that we can reproduce known results by solving the phase equation Eq. describing the massless 
charge carriers of graphene, for the case of zero-energy states formed in a Lorentzian potential V (r) = — Vb /(l + (r/d)^) 
, an analytically solved problem.l^In this case the condition for bound states when m > 0 is Vod = 2(N + m), where 
N is a positive integer. Thus, when solving Eq. 0, the phase equation for massless Dirac particles, when fc —>■ 0 
we should see the threshold value of Vod reached before step-like behavior as Vod is turned up and more confined 
zero-energy states appear. This is exactly what we find in Eig. (top). We also note that the sign of Vb is irrelevant 
for the creation of zero-energy confined states since a well for an electron is a hill for a hole and vice versa, thus 
Fig.g has mirror symmetry about the Vod = 0 axis. At the Dirac point, when the density of states vanishes yet the 
conductivity remains finite, these zero-energy states should be important as we shall see in the later on. 

Whilst screening only effects the strength of the Coulomb potential and not its characteristic decay,^ due to the 
quasi-relativistic nature of the carriers in graphene, a cut-off is necessary at the origin and the presence of an image 
charge in the metallic back gate will lead to a dipole-like 1/r^ decay at large distances, thus a convenient choice of 






























model potential is 


5 


Vir) = - — -. 

^ ^ 1 + (r/d )3 


(16) 


In Fig. (bottom) we investigate confined states with Eq. (16) and again see a characteristic threshold potential 
strength product spread Vod, followed by the signature staircase behavior of confined zero-energy states. Of course, 
compared to the exactly-solvable Lorentzian potential with its accidental degeneracies, the staircase does not share the 
same beautiful symmetries and the condition for full-confined states can be approximated by Vod « 2.63{N -|-0.72m— 
0.13). However we can now predict in realistic graphene flakes, where we would expect charged impurities to cause 
potentials similar to the type Eq. that there is the possibility for the appearance of fully-confined zero-energy 
states. 

Note, as we show in Appendix]^ by considering a general class of potential wells regular at the origin and decaying 
like some power law (which is faster than the Coulomb potential), the toy model interaction potential one chooses 
is not crucial for the actual existence of zero-energy bound states, but is important for knowing where one should 
search for them in terms of potential strength. Thus we have shown that it is possible to create electrostatic traps 
holding electrons (or holes) in the regime where the Fermi level is close to the Dirac point energy. Then it should 
be possible to release these trapped quasiparticles into the system by a small adjustment of the trapping potential 
strength, which will result in the sought-after on/off behavior. 

Experimentally, such states should be able to be detected by SPM experiments, as proposed in Ref. [36], where 
smoothly changing either the charge on the SPM tips or their distance above the graphene plane, and continuously 
holding the Fermi level at the Dirac point using the back-gate should be sufficient to see confined zero-modes. A 
network of sparse SPM tips, of radius Rup and separated in a square grid defined by an inter-tip separation s, all held 
at a distance above a metallic back gate and ft .2 ~ above the graphene plane gives a similarly behaved potential 


to Eq. (16), but in a more complicated form due to the method of images. 


^Qtip 


^(^) = fir)= {ix + js)‘^ + {y + ksf + {h 2 -hi)‘^) -{{x + jsf + {y + ksf + {h 2 + hi)"^) 

j,k=-ni 

(17) 


Fitting Eq. (16) to Eq. (|17[) by matching the functions at both their maximum and half-maximum values, one finds 


that the tip voltage Vup at which one would see such zero-modes is 


^tip — 


{Vod) 


No 


hvp 


f{0)d eRup’ 


subject to the constraint 2f{d) = /(O), 


(18) 


where we have seen from Fig. (bottom) the dimensionless parameter {Vod)N,m^o = 4.16,6.06,6.79... such that we 
are dealing witfP^ tip voltages of the order of tens of mV, Vup = 20mV, 29mV, 33mV and so on. The existence of 
confined states opens up the possibility of Coulomb blockade-type physics in graphene. Indeed, the quantum dots 
created with careful adjustment of the key parameters can be be seen to be ‘optimal traps’. An estimate of the 
charging energy of the optimal trap, using a simple disc capacitor model, shows a charging energy of the order of 
meV, thus such effects could be seen at room temperature. The effect arises due to the tightly confined nature of the 
wavefunction in optimal traps which leads to small capacitance and a significant charging energy, which is negligible 
for the usual deconfined states. 


We also show in Fig. that the zero-modes of the Dirac equation in the pote ntial Eq. (16), as predicted in the 
continuum model, are indeed present as shown via tight-binding calculations .l^SESl is found the wavefunctions have 
a ring-like structure, which ensures avoidance of any Klein tunneling effects, i.e the states are zero-energy vortices 
with m 7 ^ 0. It is most noticeable how by adjusting the parameter Vod one can go from a tightly confined state to a 
state with a highly spread probability density. We have also checked the zero-modes are indeed robust to the shape 
of the graphene flake. 

As mentioned previously, once the scattering phase shift is known a number of other useful physical properties can 
be quickly calculated. The partial cross-section Cm can be taken as 


kCm — sin {Sm)j 


(19) 


whilst the total scattering cross-section C, a measure of the length felt by oncoming particles, and the transport 
cross-section (t easily follow from the scattering phase shifts via 

. oo 2 °° 

C=^ X! sin^('^™)’ ^ X! sin2(<^m+i “'^™)> 

m——oo m— — (x> 


1/2 


(20) 
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FIG. 2: (Color online) Probability density plots of near zero-energy states confined within the model confining potential 
Eq. (161, as calculated via tight-binding methods. We show example critical V^d = 4.10 (left) and non-critical V^d = 2.60 
(right) states, as well as the associated cumulative probability plots (right). 


and we note for low-energy scattering we can take the s-wave approximation, i.e. only small m need to be considered 
as Eq. (201 is derived from partial-wave expansions where terms with high m are negligible when fc —>■ 0. Detailed 
derivations of Eqs. (19 20) via 2D elastic scattering theory adapted for Dirac fermions have been given by Novikov 
in Ref. [Ml. 




FIG. 3: (Color online) A plot of dimensionless partial scattering cross-section as a function of scaled energy k for massless 
Dirac particles incident on the model potential with 1/r^ decay given by Eq. (181, with example noncritical potential strength 


(left) Vod = 1.00 and critical (right) Vod = 2.27. We show results for m = 0,1, 2, corresponding to the solid line (red), dashed 
line (blue) and dotted line (green) respectively. 


The energy dependence of scattering cross-section has has previously been considered for a square welP^HH and 
we outline a thorough solution to this problem in Appendix We will now revisit the problem with a smooth, 
short-range potential given by Eq. which is relevant for gated structures or for hypercritical charges. We show in 
Fig. 1^ plots of dimensionless partial cross-section fcCm for the case of non-critical (left) and critical (right) parameters 
of the potential Eq. ( |l^ . We find, in contrast to the general case, at the critical potential strength the m = 0 non- 
square-integrable state goes not go to zero in fcCm as quickly as k tends towards zero, and thus is divergent in partial 
cross-section (m as energy tends towards zero. This is because to = 0 is a resonant state with a non-normalizable 
wavefunction (and a non-zero probability density at the origin r = 0) and so has such a particle has an enhanced 
likelihood of being scattered. Of course, this behavior remains in calculations of transport scattering cross-section via 
Eq. (|20l). 


We also note, in striking contrast to non-relativistic particles, at large energies the phase shift Sea is non-zero.l 


rfol 

























7 


Remarkably, 6ao is also angular momentum-independent and we find from Eq. 0 the explicit form 


(5oo = - 


y(r)dr, 


(21) 


such that for the considered potential Eq. (18) the dimensionless partial scattering cross-section in the large energy 
limit is given by displayed in Fig. Thus, the transport cross-section (t vanishes in this 

limit. 

The smoking gun of confined zero-modes in the laboratory could be via their contribution to resistivity, which in 
semiclassical Boltzmann theory can be expressed a^^ 


h 2n 
4e^ 7rn 


S 

e 


kCr 


( 22 ) 


where rig is the density of scatterers and the electron density Ue = k'^/ir. We show in Fig. I^the behavior of resistivity 
at small kd. The presence of confined zero-modes sees a much slower drop to zero than in the critical cases, such that 
when one factors in the energy dependence of rie oc there is a divergence in resistivity contribution at energies 
tending towards zero in the critical case only, p[h/Ae^] —>■ oo. It should be possible to utilize this consequent drastic 
suppression of mobility to create an off-state in graphene, which can be explored artihcially with a series of point 
gates held above the graphene monolayer. In the non-resonant case, when taking into account Ue oc we find as 
fc —?► 0 resistivity tends to a constant (depending on the number of scatterers Us) as expected. 



FIG. 4: (Color online) A plot of resistivity p, measured in units of /i/4e^, as a function of energy kd for massless Dirac particles 
scattering on the model potential with 1/r® decay. We show results for example uncritical cases Vod = 1.00,3.00, corresponding 
to the dashed line (red) and dotted line (green) respectively, and example critical cases Vod = 2.27,4.87,..., corresponding to 
the solid line (blue) and dash-dot line (purple) respectively. We set 2ns = tTOe to more clearly contrast the differences due to 
the effect of cross-section only. 


IV. DISCUSSION 

We have derived the phase equations for the 2-D Dirac equation using the VPM, suitable for use straightaway in 
scattering calculations concerning Dirac materials such as graphene. In doing so, we provide a numerical (experi¬ 
mental) proof of the Levinson Theorem for massless 2-D particles. These phase equations are of first-order, and so 
relatively undemanding computationally, and have solutions in terms of scattering phase shifts, thus other desired 
scattering properties readily follow. Applying the method to fully-confined states in graphene, we reproduce an exact 
result from the literature and go on to investigate a more physical potential again finding that a certain potential 
strengths and spreads zero-energy bound states are likely to form, which is most important when describing resonant 
scattering. 

We have also calculated the energy dependence of scattering cross-section, hnding a major distinction for the m = 0 
mode for critical potential strengths (those able to support truly bound states). In this special case, we predict a 
dominant peak in scattering cross-section at zero-energy, suggesting a high probability of being scattered, which can 
be explored experimentally via scanning probe microscopy. In an experiment with a series of point gates above the 





graphene layer, one may be able to use artificial resonant scattering to switch off the chiral tunneling effect found in 
monolayer graphene by greatly reducing the mobility of carriers. 

It should also be mentioned that the existing experiments on atomic clusters on a graphene surface can arguably be 
explained by the considered zero-energy states without involving the atomic collapse picture. These experiments deal 
with back-gated graphene which precludes the use of a long-range Coulomb potential. The measured wavefunction 
density has a pronounced ring-like structure on a scale of many lattice constants, instead of a sharp peak at the 
atomic scale as expected from ‘fall-into-the-center’ physics. Moreover, the observed features in the density of states 
are very close to zero energy (in fact, on both sides of it) instead of being a few electron volts away from this energy 
as expected for a collapsed state. 
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Appendix A: Bound states in potential wells decaying as a power law 

Let us examine the condition on bound states to form in regularized potential wells decaying asymptotically as a 
power law governed by a parameter p > 1, namely 


V{r) = -Eo 



if r < i? 
if r > i? 


(Al) 


such that the wavefunction components inside the well are simply Bessel functions and outside the well are Bessel 
functions, in a new variable ^ j product the function . Bound zero-energy states can never arise 

for TO = 0, —1, due to the requirement of some rotation to combat the Klein tunneling phenomenon.!^ Furthermore, 
some other negative angular momentum states rrie < to < — 1 are excluded depending on the power of the decay, 
for p > —2TOe. Thus, as one considers increasingly short-range interactions, one will find a wider band of ‘missing’ 
nonpositive integer angular momentum states cumulating in rrie- The eigenenergies corresponding to the square- 
integrable solutions can be found from the determinant of the matrix 


J\m\{yQR) 




= 0 , 


a = 


TO -I- p/2 

p- 1 


(A2) 


which can be solved by the usual root-finding methods for the values of VqR which are able to support confined modes. 
In the limit of large VqR >> 1, Eq. (|A2|) can be rewritten in terms of elementary functions 


cos {VqR 
cos(^ 


||TO|-f) COS {VqR - ^\m + l\ - j) 

f|a-l|-f) cos(^-flal-f) 


= 0 , 


(A3) 


and the most symmetric case p = 2 is equivalent to VqR = ^ (|to| -I- |to -I- 1| -|- 2n -|- 1), where n is an integer. 

Whilst this somewhat kinky class of potential is not as realistic as the smooth potentials considered in the main 
body of this work, it quantitatively demonstrates how bound zero-energy states are realizable in many situations for 
fast-decaying potentials. As expected, a separate analysis for the regularized Coulomb potential (p = 1) yields no 
possibility for bound zero-energy states. This conclusion is formed as the associated transcendental formed can never 
have a solution. 


Appendix B: Scattering by a circular finite potential well 


Let us consider the example of a circular finite potential welP^ given by V{r) = —VoQ{d — r), where Q{z) is the 
Heaviside step function, which allows us to obtain the following solutions of Eq. Q: XA{r) = AjnJm{{k + Vo)r) inside 
the well and XA{r) = [Jm{kr) cos(^m) — Nm{kr) sin((5m)] outside the well. Upon finding xsir) from Eq. (3b) and 
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matching both wavefunction components at the boundary r = d, we obtain the following expression for the tangent 
of the phase shift 


tan(5„^) = 


Jm{{k + Vo)(i) Jm+l{kd) — Jm+l{(k + Vo)(i) Jmjkd) 
J„((fc + Vo)d) iV„+i(fcd) - Jm+i{{k + Vo)d) N^{kd) ■ 


(Bl) 


It ca n be seen that resonances in scattering cross-section do exist at certain energies, as can be seen by substituting 
Eq. (Bl) into Eq. (19), or into the more convenient form = 1/(1 + tan((5m)“^). At energies tending to the Dirac 
point energy kC,m and, more importantly, the true partial cross-section C,m goes to zero. At large energies we find the 
expected constant value, given by kC,rn —>■ sin^(Voa)) as follows from Eq. Both of these features can be seen in 

Fig. [5] (left). 

One can find the most prominent resonances of partial cross-section for this problem, defined by k(^m = 1, which 
can only occur for some specific potential strengths, which satisfy the condition Jm{(k + VQ)d)Nm+i(kd) = Jm+i{{k + 
Vo)d)Nm{kd). If there is a solution to this equation, it is straight forward to find the resonance energy, which is the 
global maximum solution of = 0. 

In the s-wave approximation, one can see by using Eq. (Bl) in conjunction with Eq. (20) the contribution to 
resistivity, calculated via Eq. (22), is given bjP^ 


4e2 


(B2) 


which, of course, can be neglected when the concentration of scatterers and/or the scattering radius d is small. 
However, this analysis neglects the possible existence of confined zero-modes, which can be determined by the condition 
Jm{Vod) = 0. The effect of the extended m = 0 state (which occurs at the zeroes of the zeroth Bessel function, 
Vod = 2.40,5.52,...) is to see the quantity tan((5m) approach zero slowly (logarithmically) as we approach the Dirac 
point energy, such that k(^o effectively tends to a constant, producing a ‘super-resonance’ at zero-energy in pure partial 
cross-section as displayed in Fig. (right). This leads to a noticeable contribution to resistivitjESI 


h Us 1 

^ 4e2neln^(fcd) 


(B3) 


which is clearly more dominant than that found in Eq. (B2) in the s-wave approximation. 



FIG. 5: (Color online) A plot of dimensionless partial scattering cross-section fccjm as a function of scaled energy k for massless 
Dirac particles incident on the circular finite potential well, with noncritical potential strength (left) Vod = 1 and critical (right) 
Vod = 2.40. We show results for m = 0,1,2, corresponding to the solid line (red), dashed line (blue) and dotted line (green) 
respectively. 
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